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Abstract 

This paper presents the potential measurement at 1.4 TeV CLIC of the cross-section (times 
branching ratio) of the Higgs production via WW fusion with the Higgs subsequently de¬ 
caying in ZZ*, <j(Hv e v e ) x BR(H — > ZZ*), and of the Higgs production via ZZ fusion with 
the Higgs subsequently decaying in bb, a (He + e~) x BR(H -> bb). For the H -> ZZ* de¬ 
cay the hadronic final state, ZZ* —> qqqq, and the semi-leptonic final state, ZZ* — > qql + l~, 
are considered. The results show that a(Hv e v e ) x BRiH —> ZZ") can be measured with 
a precision of 18.3% and 6% for the hadronic and semi-leptonic channel, respectively. 
a(H e + e ) x BR(H —> bb) can be measured with a precision of 1.7%. This measurement 
also contributes to the determination of the Higgs coupling to the Z boson, gH zz - 
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3 The CLIC ILD detector model 


1 Introduction 

The Compact Linear Collider (CLIC) is a proposed high-luminosity linear e + e~ collider planned to be 
implemented in stages with centre-of-mass energy, y/s, of 350 GeV (or more), 1.4 TeV and 3 TeV. One 
of the main aims of CLIC would be the high precision measurement of the Higgs boson properties [1,2]. 

This paper focuses on the measurement contributing to the determination of the Higgs coupling to the 
Z boson, gn zz , at 1.4 TeV CLIC. 

At s/s = 1.4 TeV, the dominant Higgs production process is the WW fusion, with ~ 370000 expected 
events in 1.5 ah 1 of data. This would lead to the measurement of the relative coupling of the Higgs boson 
to the W and Z bosons at the percent level, providing a strong test of the Standard Model prediction for 
8H ww /gH zz = cos 2 0w, where Qw is the Weinberg angle. 

The subleading Higgs production process is the ZZ fusion that, with ~ 10% the cross section of the 
WW fusion, would give ~ 37000 expected events in 1.5 ab -1 of data and it would provide access to 
complementary Higgs boson coupling, gH zz , at a percent level. 

In this paper, we discuss the measurement of the o(Hv e v e ) x BR(H —> ZZ*) where Higgs is produced 
via WW fusion, and the fully hadronic ZZ* — > qqqq, and semi-leptonic ZZ* — > qql l final states are 
considered. For the Higgs production in ZZ fusion, the H —>• bb decay is analysed. 

Both analyses are performed using the CLIC_ILD detector concept [3], assuming a total integrated 
luminosity of 1.5 ah 1 and unpolarised beams. 

2 Simulation and Reconstruction 

The Higgs production through WW and ZZ fusion is generated in Whizard 1.95 [4] , where a Higgs mass 
of 126 GeV is assumed. The background events are also generated in Whizard. using Pythia 6.4 [5] to 
simulate the hadronisation and fragmentation processes. The CLIC luminosity spectrum and the beam- 
induced processes are obtained from GuineaPig 1.4.4 [6]. All events are simulated with unpolarised 
beams. 

The interactions with the detector are simulated using the CLIC_ILD detector model within the Mokka 
simulation package [7], based on GEANT4 [8]. 

Events are reconstructed using the particle-flow technique, implemented in the Pandora algorithm 
(PFA) [9, 10]. The k t algorithm [11], implemented in FastJet [12], is used in the exclusive mode to cluster 
the jets of each event. The LCFIPlus package [13] is used for the identification of charm and beauty jets. 
The IsolatedLeptonFinder Marlin processor [14] is used to identify leptons. The TMVA package [15] is 
used for the multivariate classification of signal and background events using their kinematic properties. 
The simulation, reconstruction and analysis are done with ILCDIRAC [16], 

3 The CLICJLD detector model 

The CLICJLD detector is based on the ILD detector concept [17] for ILC, modified according to spe¬ 
cific experimental conditions at CLIC [1]. The main tracking device of CLICJLD is a Time Projection 
Chamber providing a point resolution in the r(j) plane better than 100 pm. The precision physics at CLIC 
require a vertex-detector system with excellent flavour-tagging capabilities through the measurement of 
displaced vertices. The vertex detector is based on ultra-thin hybrid pixel sensors technology, and uses 
power-pulsing and air-flow cooling to minimise the material budget. The CLICJLD detector concept 
is based on fine-grained electromagnetic and hadronic calorimeters (ECAL and HCAL), optimised for 
particle-flow techniques. Both calorimeters are within a 4 T solenoidal magnetic field. 
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4 Measurement of <j(Hvv) x BR{H —r ZZ*) 


V e 



Figure 1: Feynman diagram of the Higgs production in WW fusion and the subsequent Higgs boson 
decay to a pair of Z bosons. 

4 Measurement of a(Hvv) x BR{H -* ZZ*) 

4.1 Event samples 

The WW fusion process has the largest cross section for Fliggs production at 1.4 TeV CLIC. The cross 
section for e + e —r Hvv is 244 fb. The branching fraction for the H —> ZZ* decay is 2.89% [18]. The 
Feynman diagram for the process e e" — y Hvv — > ZZ*vv is shown in Figure 1. 

In this analysis, two final states of the ZZ* decays are studied: the fully hadronic final state, ZZ* — > 
qqqq, with a branching ratio of 49%, resulting in an effective cross section of 3.45 fb and the semi- 
leptonic final state, ZZ* — > qql ' / , with a branching ratio of 10% and an effective cross-section of 0.72 
fb. 

In Table 1, the full list of the signal and background processes is given with the corresponding cross- 
sections. 

The main background, characterised by the same final state as the fully hadronic signal final state, is 
given by the e~e + Hv e v e —> WW v ( ,v,, —? qqqqv e V e process. Other important backround processes 
due to their large cross section are e + y^t qqqqv, 77 —> qqqq and t ,+ 7 ^ qqqqe. They can be substantially 
reduced by requiring high-p T jets. Other minor background processes can be discriminated from signal 
events using a MVA event classifier. 

4.2 Analysis strategy 

The following section describes the physics object identification and the selection requirements applied 
in the analysis. 

For the semi-leptonic final state, the first step of the physics object identification is searching for 
isolated leptons (electrons or muons). Exactly two leptons are required, otherwise the event is rejected. 
Then, all particles in the event not identified as leptons are clustered by the k, algorithm into two jets 
with R = 1.0. 

Instead, for the hadronic final state, the event is directly clustered by the k, algorithm into four jets 
with R = 1.0. 

Next, for both final states, flavour-tagging is performed and a preselection based on kinematics vari¬ 
ables is applied. Finally, a MVA event selection based on the BDT classifier is performed to obtain the 
final results. 

4.2.1 Lepton identification 

The first step of the semi-leptonic analysis is to identify and reconstruct leptons from ZZ* decays. This 
section describes how leptons are distinguished from all other reconstructed particles in the event. Isol¬ 
ated leptons are identified using a combination of track energy, calorimeter energy and impact parameter. 

Muons and electrons are required to have a track energy of at least 7 GeV. The impact parameter of a 
track describes the perpendicular distance between the track and primary vertex (PV), at the track’s point 
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Table 1: List of considered processes with the corresponding cross-sections. Background processes 
marked by f, where jets and leptons arc not produced by a Higgs boson, arc generated with a 
cut on the Higgs mass set to m#=12 TeV. Because of the large cross section of the processes 
marked by *, generator level cuts are introduced to limit the number of events that needed to be 
simulated and reconstructed. The invariant mass of the two jets, m qq , is required to be greater 
than 50 GeV and invariant mass of four jets, m qqqq is also set to be greater than 50 GeV. The 
cross sections for all processes with photons in the initial state include both processes with the 
Beamstrahlung and with Effective Photon Approximation (EPA) photons. Cross sections for 
the processes e + y —> qqv, e + y -A qqqqV, e + y —» qqqqe, and e + y -A qqe represent the sum of 
processes with the electron and the positron in the initial state. 


Signal process 


o(Jb) 

e~e+ -A- Hvv,H -A ZZ\ZZ* - 

-7- qqqq 

3.45 

e~e + -A Hvv,H -A ZZ\ZZ* - 

A qql + l~ 

0.72 

Common background 


o{fb) 


e e + -A qqv e V e 


788 

e~e + -A qqqqV e V e 


24.7 

e~e + -A Hvv,H — 

► WW, WW -A qqqq 

27.6 

e~e + -A Hvv,H — 

> bb 

136.94 

e~e + -A- qq 


4009.5 

e~e + -A- qqqq 


1245.1^ 

e~e + -A qqqql + l~ 


71.7 

e~e + -A qqqqlv 


115.3 

e~e + -A Hvv,H — 

yzz*,zz* -Ai+ri + i~ 

0.08 


-a Hvv,H -A ZZfZZ* - 

A qqqq specific background 

o(fb) 

e ± y -A qqqqV 


338.5 t 

77 -> qqqq 


30212*f 

e + 7 ^ qqqqe 


2891 

e~e + -A Hvv,H -A ZZ*,ZZ* - 

qql + l~ 

0.72 

e~e+ -A Hvv,H -A ZZ\ZZ* - 

A qql + l~ specific background 

<y(fb) 


e e 
e ± y 
e ± y 

77 - 

e~e^ 


—> qql + l~ 

A qqv 
A qqe 
qq 

—> Hvv,H -A ZZ*,ZZ* -A qqqq 


2725.8' 
37125.3 *f 
63838.8*f 
112038.6 *f 
3.45 


of closest approach to the PV. It can be decomposed into longitudinal (Zo) and radial (do) components, 
which combine to give the impact parameter in 3 dimension (Ro) 

Ro = \J Zq 2 + do 2 • ( 1 ) 

In the analysis, an impact parameter smaller than 0.02 mm is required. 

The ratio (RCAL) of the energy deposits in the ECAL and in the HCAL 


RCAL = Eecal/ ( Eecal + Ehcal ) (2) 

is a powerful discriminating variable. 

Since electrons are mostly contained within the ECAL, they have a peak at RCAL=1. Muons deposit 
a minimal amount of ionisation energy throughout the calorimeters, and have a peak at RCAL=0.2. 
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Therefore, in order to reduce the miss-identification lepton contribution, RCAL is required to be either 
larger than 0.9 or in the range 0.025-0.300. 

Applying these selection criteria, 74% of electron and muon pairs from Z decays are correctly identi¬ 
fied. 

4.2.2 Preselection 

Leptons and jets are paired to give the Z bosons contribution. Since mu < 2 mz, one pair is required to 
have a mass consistent with mz (on-shell Z boson), while the second pair is required to form the off-shell 
Z boson. The preselection cuts require: 

• on-shell Z boson: 45 GeV < mz <110 GeV, 

• off-shell Z boson: mz* < 65 GeV, 

• Higgs invariant mass: 90 GeV < mu < 165 GeV, 

• the distance value between the two closest jets (jet transition cuts): — logy 34 < 3.5, — logy 23 < 3.0, 

• visible energy: 100 GeV < E VIS < 600 GeV, 

• missing transverse momentum p™ ss > 80 GeV, 

• In order to reject H —>• bb decays, the event is forced into a two-jet topology and the flavour-tag is 
applied to the two jets. Events where one or both jets have a b-tag probability, P(b)^ et \ greater 
than 0.95 are rejected. 

The jet transition values, — log 3:34 and — log v’ 23 , are used in the preselection to discriminate signal 
from lower jet multiplicity backgrounds. 

In addition, for the semi-leptonic final state exactly two isolated leptons are required. 

The signal preselection efficiency is 30.2% for the fully-hadronic final state and 74% for the semi- 
leptonic final state. The preselection efficiency in the fully-hadronic final state is relatively low due to 
p™ ss > 80 GeV cut in order to suppress the e + y qqqqv and e + y —t qqqqe processes. 

4.2.3 MVA event selection 

After the preselection, a MVA event selection based on the BDT classifier is applied in the analysis. 

For the fully hadronic final state, the following 11 sensitive observables are used for the classification 
of the events: mz, mz*, mu, — logy 34 , — logV23, Evis, Pt, P{b)i etl , P(b)' et2 , P(c)-i etl and P(c)- 7 ^ 2 ; while 
for the semi-leptonic final state the following 17 observables are used: mz, mz*, m/+/-, m qq , mu, E v i s , 
-logy 34 , -logy 2 3 , -logyn, P(b) jet p P(b) jet \ P{c) iet \ P(c) jet2 , p T , 0 H , E vis - E„ and N PFO - These 
variables are defined as follows: Nppo is the number of all particle-flow objects in one event, m/+/- is the 
invariant mass of the two selected leptons, m qq is the invariant mass of the two selected jets, mz is the 
invariant mass of the on-shell Z boson, mz* is the the invariant mass of the off-shell Z boson, mu is the 
invariant mass of the Higgs candidate, E v - ls is the visible energy of the event, E v - ls — Eu is the difference 
between the visible energy in the event and the Higgs visible energy, — logy 34 ,— logy23 and — log)’12 are 
transition variables, P(b)E tl and P(b) ' eh - are b-tag probability of the jets, P(c)i et ' and P(c ) /> ' 2 are c-tag 
probability of the jets, pr is missing transverse momentum and % is polar angle of the Higgs candidate. 

The BDT is trained on H —> bb, e e + —> qqv e V e and e ± y qqqqv background samples and the signal 
sample for the fully hadronic final state, while it is trained on all the background sample and signal for 
the semi-leptonic final state. 

In both final states, the BDT cut maximising the significance is chosen, giving an overall efficiency of 
18% and 33%, for the fully-hadronic and semi-leptonic final states, respectively. Figure 2 (left) includes 


5 
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■ Signal 3 

■ H_others 8 

■ qqqqlnu 

■ qqnunu 

■ qqqqnunu 

■ Hbb 

■ HWWqqqq 
qqqqll 

■ qqqqnu 



Figure 2: The Higgs invariant mass distributions after preselection (left) and after MVA selection (right) 
for the fully-hadronic final state. 
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Figure 3: The Fliggs invariant mass distributions after preselection (left) and after MVA selection (right) 
for the semi-leptonic final state. 


all events that pass the preselection, while Figure 2 (right) shows all events passing the BDT selection 
for the fully-hadronic final state. The same is shown for the semi-leptonic final state in Figure 3. All 
samples are normalised to the integrated luminosity of 1.5 ah 1 . 

The final results and the statistical uncertainty for the a(Hv e v e ) xBR(H —>ZZ*) measurement are 
reported in Section 6. 

5 Higgs production in ZZ fusion and a(He + e~) x BR(H —> bb) 
measurement 

5.1 Event samples 

The ZZ fusion process has a 10 times smaller cross section for Fliggs production at 1.4 TeV CLIC than 
the WW fusion process. The characteristic signature of the ZZ fusion process is two scattered beam 
electrons with a large pseudorapidity separation, plus the Higgs boson decay products. The cross section 
for e + e~ -» He ~e + is 25 fb. In this analysis, the scattered beam electrons are required to be fully 
reconstructed, and the final state H — > bb is considered. The branching fraction for the H -x bb decay is 
56.1%. The corresponding cross-section for the signal is 13.74 fb. 

In Figure 4 the Feynman diagram for e + e~ —> He e 1 —> bbe~e 1 is shown. In Table 2 a full list of the 
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5 Higgs production in ZZ fusion and o(He l e ) X BR(H —> bb) measurement 



Figure 4: Feynman diagram of the Higgs production in ZZ fusion and the subsequent Higgs boson decay 
to a pair of b quarks. 


Table 2: List of considered processes with the corresponding cross-sections. 

Process 

e~e + —> He~e + ,H —> bb 
e~e + —> qql + l~ 
e~e + —> qqqq 
e~e + —> qqqql + l~ 
e~e + -* tt 
77 qqe~e + 


<y(fb) 

13.74 

2727 

1328 

71.7 

135.8 

6 


signal and background processes is given with the corresponding cross-sections. 

The main backgound is e~e + —> qql + l~. It can be reduced requiring b-tagged jets. Other background 
processes give very small contributions after the preselection. 

5.2 Analysis strategy 

This section describes the analysis method for the reconstruction of Higgs candidates produced in ZZ fu¬ 
sion undergoing subsequent decay into bb. Events are clustered into a 4-jet topology using a k t exclusive 
clustering algorithm with R = 1.0. For a well-reconstructed signal event, two of the resulting ‘jets’are 
expected to be reconstructed electrons, and the remaining two jets are from the Higgs decay to bb. The 
discrimination between signal and background events is based on pre-selection cuts and a multivariate 
likelihood analysis. 

5.2.1 Preselection 

The preselection requires two oppositely-charged electron candidates, separated by | Ar/1>1, each with E 
> 100 GeV. A further requirement is that at least one of the two jets associated with the Higgs decay has 
a b-tag value greater than 0.65. These cuts lead to signal efficiency of 19.3% resulting in an effective 
cross-section of 2.6 fb. The major acceptance loss comes from the geometrical effect of electrons falling 
outside the detector (Figure 5). 

The main background process after the preselection is e~e + —> qql l selected with a 0.2% efficiency 
resulting in an effective cross-section of 6.44 fb. 

5.2.2 MVA 

As a next step in the analysis method, a classification and selection based on a multivariate data analysis 
is performed. A relative likelihood classifier is constructed using four variables that provide separation 
between signal and background: 

• opening distance between the reconstructed electrons: AR, 
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6 Results 



Figure 5: Scattered beam electrons at the generator level for different centre-of-mass energies (red ar¬ 
rows indicate detector acceptance). 



Signal likelihood 


Figure 6: Likelihood distribution for signal (red) and background (blue). 


• recoil mass of the event determined from the momenta of the reconstructed electrons: m rec , 

• jet transition variable: — logj 34 , 

• the invariant mass of the two jets associated with the Fliggs decay. 

The resulting likelihood distribution is shown in Figure 6 and gives good separation between signal 
and background. The result for the statistical uncertainty is reported in Table 3. 


6 Results 

The results of the measurements of o(He + e~ ) x BR(H — > bb), g(H v t ,v e ) x BR(H —> ZZ* H> qqqq) and 
o(Hv e v e ) x BR{H —> ZZ* —> qql + l ) are shown in Table 3. All measurements are simulated at 1.4 TeV 
CLIC collider with unpolarised beams. 

The relative statistical uncertainty is 18.3% and 6% for the hadronic and semi-leptonic ZZ* decays, 
respectively. They are dominated by the limited signal statistics and the presence of large backgrounds 
in the ZZ* —> qqqq measurement. 

The obtained results are included in the the global fit to contribute to the Higgs to Z coupling, g// Z2 , 
and to the total Higgs width F//. 
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Table 3: Summary of the results of the analysis of the o(Hvv) x BR{H —>• ZZ*) and o(He + e~) x 
BR(H —> bb) measurements at 1.4 TeV CLIC with unpolarised beams. N s is the number of 
final signal events, £ s represent overall signal efficiency of the measurements and 5 (a x HR) is 
statistical uncertainty on the cross section times branching ratio of the measurements. 



Hw -> ZZ'.ZZ' —>■ qqqq 

Hvv -> ZZ*,ZZ* -> qql+l~ 

He + e —> bb 

N s 

931 

297 

748 

£ s 

18% 

33% 

19.3% 

5(a x BR) 

18.3% 

6% 

1.7% 


These results may be improved including tau leptons in the semi-leptonic analysis and trying to further 
improve the lepton pair efficiency. 

A further improvement may come from the beam polarisation: if 80% left-handed polarisation of the 
electron beam is assumed during the entire operation time at 1.4 TeV, the WW fusion Higgs production 
cross-section would be enhanced by a factor 1.8. 

The statistical uncertainty of the measurement a(He + e~ ) x BR(H —> bb) is 1.7%. This measurement 
is proportional to g 2 Hzz , ■ gjj _/T // and the result is included in global Higgs fit to contribute to the gn z/ 
determination. 
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